Brownian motion in superfluid "^He 



O 
(N 



Xiao Li^,* Ran Cheng^, Tongcang Li, and Qian Niu^ 
Department of Physics, The University of Texas at Austin, Texas 78712, USA 

We propose to study the Brownian motion of a classical microsphere submerged in superfluid ^He using the 
recent laser technology as a direct investigation of the thermal fluctuations of quasiparticles in the quantum fluid. 
By calculating the temperature dependence of both the friction coefficient and the strength of the random force, 
we show that the resonant mode of the fluctuational motion can be fully resolved by the present technology. 
Contrary to the previous work, it is found that the roton contribution is not negligible, and it even becomes 
dominant when the temperature is above 0.76 K. 
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Since its discovery in 1827, Brownian motion has at- 
tracted broad interest and become a major thrust of statistical 
physics'"^. Driven by thermal fluctuations among the con- 
stituent particles, the Brownian motion gave the first direct 
evidence of the particle nature of the background fluid. How- 
ever, at extremely low temperatures, quantum effects of the 
fluid particles become important and certain kinds of quantum 
fluid emerge, where the system can be described by a macro- 
scopic wave function with long-range phase coherence. In 
such a case, it is the elementary excitations that play the role 
of the fluid particles. In thermal equilibrium, those excitations 
form a gas of quasiparticles that behaves in a similar way to 
the ordinary atomic gas. They are able to create pressure and 
resistance when colliding at the interface with a solid. How- 
ever, their unusual dispersion relations are quite different from 
those of ordinary particles. A good example is the superfluid 
^He where two kinds of particle-like excitations, phonons and 
rotons (see Fig. 1), comprise the viscous part of the fluid^**. 
A natural question arisen is whether these quasiparticles can 
create Brownian motion. To date, this question has not been 
addressed by any experiment, and the only theoretical study 
by Balazs'^ has ignored the roton contribution without justifi- 
cation. Moreover, Balazs' prediction was far beyond the ex- 
perimental capability for two reasons: first, the amplitude of 
the Brownian particle was too small to be monitored at such 
low temperatures; second, the required observation time was 
too long due to the tiny restoring force provided by the quartz 
fiber. But without a fiber, on the experimental side, no object 
could be suspended in superfluid ''He for its small density. 

A recent progress in laser technology, however, removes 
these limitations. It not only promotes the position and time 
resolutions to a very high level'"'", but also is capable of 
trapping a classical microsphere by laser beams of high inten- 
sity that provides a sufficiently strong restoring force, roughly 
about 10^ ~' 10^ times higher than a quartz fiber. In addi- 
tion, the superfluid medium is transparent to the laser thus the 
thermal equilibrium of the fluid will not be destroyed by the 
beams. This technical development therefore makes it possi- 
ble to test the Brownian motion experimentally in superfluid 
''He, or potentially, any quantum fluid that does not interact 
with the laser. Therefore, a complete study of the Brownian 
motion in superfluid ''He, incorporating contributions from 
both the phonon and the roton, will shed light on both the- 
oretical and experimental investigations. 



In this paper, we use the Langevin equation to describe the 
dynamics of the Brownian particle with the relevant parame- 
ters taken from recent experiments'". The friction coefficient 
and the strength of the random force are evaluated by the ki- 
netic theory, and expressed as functions of the temperature. 
We conclude that the Full-Width-at-Half-Maximum (FWHM) 
and the mean square of the resonant mode of the fluctuational 
motion are able to be directly measured by the present exper- 
imental techniques'"". Interestingly, while the roton contri- 
bution to the above two quantities is negligible at low temper- 
atures, it becomes dominant above 0.76 K. 
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FIG. 1. Quasiparticle spectrum in superfluid He. We have the lin- 
ear phonon excitations at small momenta and the roton excitations 
at higher momenta. Rotons with negative and positive slopes are 
labeled by _R_ and R+ respectively. The dashed lines represent 
the analytical expressions to be used to approximate the spectrum: 



e{p) = CsP and s{p) 
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where Cs 



239 m/s. 



m = 1.06 X 10"^'^ kg, and p,„t/h = 1.92 A. 



We consider a microsphere with mass M held in super- 
fluid ^He by a laser trap with harmonic angular frequency 
cjQ- Quasiparticles excited by thermal fluctuations in the su- 
perfluid will create random forces on the microsphere, result- 
ing in the Brownian motion around its equilibrium position. 
The motion of the microsphere can be well described by the 
Langevin equation: 
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where r denotes the position of the ball, 7 represents the fric- 
tion to be estimated below, and F{t) is the random force that 



drives the Brownian motion. It is apparent that different spa- 
cial components of r are independent in Eq. (1). Thus in 
the following we shall only focus on the z-component of the 
Brownian motion. As a consequence, we now consider an ef- 
fective plate with area a = irr"^ to represent the microsphere, 
upon which the problem reduces to one dimension. 

We will solve the Langevin equation under the following 
assumption; 



< 27r/cJo < ■ 



(2) 



which says that the observation time t far exceeds the typical 
period of free oscillation, and the latter is also much larger 
than the time interval Tc between two adjacent collisions from 
the quasiparticles'. We decompose the random force into a 
Fourier sum Fz{t) = X)n(^" cos w„t + i?„ sinaj„<), where 
ujn = 27rn/r, with n taking integer values, and A„ = 
T Jo ^^ (^) '^°^ ujktdt, Bn = ^ Jq Fz (t) sin w„td<. It is worth 
mentioning that the Fourier decomposition makes sense only 
when we regard F^ (<) as a periodic function of time with pe- 
riod T. By expressing the z-component of the displacement 
as z{t) = J2n Zn{t), Eq. (1) can be solved in the frequency 
domain where 
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What subject to direct experimental verifications are the 
FWHM of the peak and the mean square amplitude of the 
resonant mode at aj„ = cjq, which equal to 7/A/ and 
(Ag + i?Q)/27^Ci;g respectively, where ( ) denotes the ensem- 
ble average. The latter one can be brought into a more use- 
ful form if we make the following considerations. Because 
the quasiparticle density is dilute, two adjacent collisions can 
be considered uncorrected, and thus we will adopt the white 
noise assumption that {F:^{ti)Fz{t2))T = 'D{T)5{ti - ^2), 
where 2?(T) is a function of the temperature alone. Also, 
by regarding r as sufficiently large (According to Ref. 10, 
Wo ~ 27r X 10'^ Hz, thus r = 1 s meets the requirement), 
we convert the summation "^f^to the integral ^ J duj. Then 
we obtain from Eq. (3) the mean square of the resonant mode. 
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Therefore, what we need to evaluate are j{T) and 2?(r) as 
functions of the temperature. Although the equipartition rela- 
tion P/47 = kBT/2^''^ obviates the need to calculate both of 
them, we still do so for strictness. 

The friction 7 originates from the imbalance of the forward 
and backward scattering when an object has an instantaneous 
velocity with respect to the background fluid. Suppose Act 
is a given area on the front side of the effective plate with 
its normal taken to be the z-axis. When the plate moves 
along z direction with velocity Vz, the number of quasipar- 
ticles within momentum interval [p, p + dp] that is able to 
collide on Act during time At (r^ <C Af <C 2tt/llIo) is given 
by lS.a ISid^ p\v z — Uz\Np/h^, where Uz ~ de{p)/dpz is the 
group velocity of the quasiparticle, and Np is the number of 



particles in each unit volume with vector momentum p. If 
each collision transfers a momentum 6pz to the plate which 
will be specified below, the resulting force from the front is 
the average of the total momentum transfer during this time 
interval divided by At, 
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where ct is the total area of the plate. Similarly, force Fi, from 
behind is given by the same expression but the range of Uz 
should be Vz < Uz < 00. Addition of F/ and Fb gives the net 
resistant force F,. = Ff + Fjy = ^Vz + 0{v1), with the linear 
term in Vz being the friction. 
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FIG. 2. Reflections and transmissions of quasiparticles on the in- 
terface separating the effective plate and the superfluid ''He. This 
diagram only show the case when a phonon with momentum p'j""^ in- 
cident on the interface with four outgoing channels: elastic reflection 
to a phonon pi with possibility Rn, to a R- roton p2 with possibil- 
ity i?i2 , and to a _R+ roton pz with possibility _Ri3 ; inelastic collision 
that creates a phonon p4 into the solid with possibility i?i4. Arrows 
on each line represents the group velocity of that quasiparticle. Note 
especially that the group velocity for R- is in the opposite direction 
of its momentum. 

Specific computation of 7 along this line needs detailed in- 
formation on 5pz. Let us now focus on an individual colUsion 
that respects the momentum and energy conservations, 

Pz + Mvz = {pz-5pz) + Mv',, 



e{p) + -Mvl = e(vV 



2pz5pz + {5pzY) + -Mv 
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and we regard the mass of the Brownian particle M as suffi- 
ciently large (A/ = 2.8 x 10^^** kg experimentally) so that 
each process is a hard wall collision. Problem arises when the 
energy of the incoming quasiparticle e(p) lies in the region 
(A^,,^^^^, A^^^.^^J, for which the above equations admit three in- 
equivalent solutions of 5pz- That is to say, for example, for 
a phonon carrying definite momentum and energy incident on 
the plate, it has three different outgoing channels labeled by 
i = 1, 2, 3 representing the phonon, i?_ and i?_|_ rotons re- 
spectively. We may also add the possibility of inelastic col- 
lision that creates phonons in the Brownian particle labeled 



as i = 4 '^•'^. This multi-channel collision process is illus- 
trated in Fig. 2. The momentum transfer 5pz then depends 
on the transition probability Rij connecting the i-th and the 
j-th channels. However, it is remarkable that the final expres- 
sion of 7 turns out to be independent of Rij as if there were 
no inter-channel transitions. We omit the detailed argument of 
this result, as it shares similar logic with the problem of quasi- 
particle pressure in superfluid ''He'-, where all inter-channel 
transitions mutually cancel. After some manipulations, both 
the phonon and the roton contributions can be brought into the 
same form as 
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Insertion of the dispersion relations e{p) = CsP and e{p) = 

A + 2^ (S^^ FiS- 1) yields the contributions from the 
phonon and the roton respectively. 
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from the fluctuation of the random force exerted on the 
Brownian particle. Mathematically, the fluctuation is em- 
bodied in the statistical deviation of the quasiparticles dis- 
tribution n{r,p) in the six-dimensional phase space which 
satisfies Np = J cfirn{r,p). Neglecting the influence on 
n{r,p) of the scattering among quasiparticles, we again 
assume that n{r,p) on different phase points are inde- 
pendent so that {n{r,p)n{r',p')) ^ (5''(r — r')6^{p — 
p'), i.e., only contributions from the same phase point are 
kept. Then with the help of the Bose relation {{Np — 

Np)^) = Np + Np , we know after some manipulations 
the total fluctuation of the momentum transfer during Ai is: 
(G^L, = aAt J d^p/h^\u,\{Sp,)^{Np + Np"^). The corre- 
sponding fluctuation of the random force should be {F'})^^ = 
{Gl)^J{At)'^ and the quantity of central interest is T> = 
(F^ )^j Ai, which gets rid of the Ai dependence. 
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Again, by inserting the dispersion relations, a straightforward 
calculation leads us to: 



and the total friction coefficient is 7(T) — jph{T) + 7rot(7')- 
For a Brownian particle with M = 2.8 x 10^'^ kg'*', the 
FWHM ^{T)/M is depicted in Fig. 3. The crossing point of 
7ph(T) and 7i-ot(r) is at 0.76 K. While the roton contribution 
is negligible at low temperatures, it becomes dominate above 
0.76 K. We see from the figure that the typical FWHM ranges 
from 10 Hz to 1 kHz, and that the departure from pure phonon 
contribution above 0.76 K is of the order of 1 kHz. These are 
well within the capability of current experiments where the 
resolution of frequency is down to 1 Hz. 
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FIG. 3. The temperature dependence of the FWHM of the mean 
square displacement spectra. The thick black line denotes the total 
width, while the dashed and dash-dotted lines are contributions from 
the phonon and the roton respectively. The effective area of the plate 
is taken to be 97r/4 /im^. 

Having fully evaluated the friction coefficient 7(T), now 
we turn to the more involved quantity P(T) which comes 
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and V{T) = V^T) + V,ov{T). We mention in passing that 
the equipartition relation 2?/47 = kBT/2 does hold sepa- 
rately for phonon and roton excitations. 

Equipped with the friction coefficient and the fluctuation of 
random force, we are able to evaluate the temperature depen- 
dence of the resonant mode, which serves as another quantity 
for direct experimental test. In view of the equipartition rela- 
tion, the resonant mode in Eq. (4) is 



2kBT/Lo'^7r 
7ph(T)+7.o,(r)' 



(12) 



This is plotted in Fig. 4, where 0.76 K is again identified as the 
turning temperature. The higher the temperature, the lower 
the amplitude of the resonant mode is. Fortunately, the lowest 

-3 



\/{z'j^) seen from the plot is about 10 -^ ~^ 10 ■^ nm/-\/Hz, far 
beyond the experimental resolution 3.9 x 10^^ nm/\/Hz'". At 
low temperatures, however, the resonant mode is dominated 
by the phonon excitation and diverges as T^'^. This seem- 
ingly counterintuitive result is resolved when we consider the 
assumption made in Eq. (2), which has been discussed in 
Ref. 9. As the temperature goes down, Tc sharply increases 
and the above assumption becomes invalid, where a new the- 
ory is required. Fortunately, for a typical angular frequency 
ujo ^ 2tt X 6000 Hz, Eq. (2) holds until the temperature is 
lowered to lO"*" K. 

Before conclusion, two further comments are in order 
First, the motion of the Brownian particle will transfer ki- 
netic energy to the background fluid. As a consequence the 
mass M in Eq. (1) should be the effective mass instead of the 
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FIG. 4. The temperature dependence of the square amphtude of the 
resonant mode. The thick hne represents the result of Eq. (12). In 
comparison, we also plotted the result in existing literatures by the 
dashed line, where the roton contribution is absent. The dash-dot 
line depicts the case if the roton contributes alone. Parameters relat- 
ing to experiment are the same as those in Fig. 2, and the resonant 
frequency uq is taken to be 27r x 6000 Hz. 



bare mass toq. Below T =1 K, it can be approximated as'"*: 

1 4 ^ 

M — mo + —ps X oTT?"^, where ps is the superfluid density 

and r is the radius of the microsphere. For the microsphere 
with a diameter of 3 /xm and mass of mo ~ 2.8 x 10^'^ kg. 



the correction term is roughly about 5%. 

Moreover, the white noise assumption on the random force 
as well as the delta correlated distribution n{r,p) implies the 
independence of different quasiparticle modes. This is quite 
reasonable when the density of quasiparticles is dilute, as the 
relative strength of quasiparticle scattering is proportional to 
its square. At roughly T =1 K, the percentage of the normal 
fluid formed by the quasiparticles is less than 5%'^, and we 
would expect a negligible effects of quasiparticle scattering. 

In conclusion, we have studied the Brownian motion of 
a classical microsphere driven by thermally excited quasi- 
particles in superfluid "^He. Contrary to previous work, we 
claim the importance of contributions from both the phonon 
and the roton excitations, and find the turning temperature of 
their relative importance at 0.76 K. More importantly, the two 
predictions we give on the FWHM and the resonant mode are 
able to be tested in current experiments. Generalization to 
other types of quantum fluids are left for future inquires. 
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